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Abstract 

The EKG or electrocardiogram sequence is defined by a(l) = 1, a(2) = 2 and, for n > 3, a(n) is 
the smallest natural number not already in the sequence with the property that 
gcd{a(ra — l),a(n)} > 1. In spite of its erratic local behavior, which when plotted resembles 
an electrocardiogram, its global behavior appears quite regular. We conjecture that almost all a(n) 
satisfy the asymptotic formula a{n) = n(l + 1/(3 log n)) + o(n/ log n) as n — > oo; and that the 
exceptional values a{n) = p and a(n) = 3p, for p a prime, produce the spikes in the EKG sequence. 
We prove that {a(n) : n > 1} is a permutation of the natural numbers and that c\n < a{n) < C2n 
for constants c\, ci- There remains a large gap between what is conjectured and what is proved. 

1. Introduction 

Consider the sequence defined by a(l) = 1, a(2) = 2 and, for n > 3, a{n) is the smallest natural 
number not in {a{k) : 1 < k < n — 1} with the property that gcd{o(n— l),a(n)} > 2. This sequence 
might be called a greedy gcd sequence, but because of its striking appearance when plotted we 
will name it the EKG (or electrocardiogram) sequence — see Figures |l|, 0. It was apparently first 
discovered by Jonathan Ayres [Ayres 2001] and appears as sequence [A064413 in [Sloane 2001]. The 
first 30 terms are 

1 24639 12 8 10 5 
15 18 14 7 21 24 16 20 22 11 
33 27 30 25 35 28 26 13 39 36 ••• 

Although the local behavior is erratic, plots of the first 1000 or 10000 terms show considerable 
regularity (see Figures |, | in Section 4). 

The EKG sequence has a simple recursive definition, yet seems surprisingly difficult to analyze. 
Its definition combines both additive and multiplicative aspects of the integers, and the greedy 
property of its definition produces a complicated dependence on the earlier terms of the sequence. 
Indeed, it is not immediately obvious whether it contains all positive integers, but we show this is 
the case — the EKG sequence is a permutation of the positive integers. 
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ure 1: Plot of a(l) to a(100), with successive points joined by lines. 
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igure 2: Terms 800 to 1000, with successive points joined by lines. 



In comparing Figures 

and |, one is reminded of the contrast between the irregular plot of tt(x) 
(the number of primes < x) for x < 100 and the very smooth plot for x < 50000 as shown in Don 
Zagier's lecture on "The first 50 million prime numbers" [Zagier 1977]. This is not a coincidence, 
as we will see, because (experimentally) the spikes in the EKG sequence are associated with the 
primes arranged in increasing order. However, the spacings between spikes are not the same as the 
spacings between consecutive primes. 

Although the EKG sequence itself seems only to have been proposed recently, in the early 
1980's Erdos, Freud and Hegyvari [Erdos et al. 1983] studied properties of integer permutations 
with restrictions placed on allowed values of greatest common divisors of consecutive terms. 

In Section 2 we derive a number of basic properties of the EKG sequence, and prove that it is a 
permutation of the natural numbers N. An efficient algorithm for computing the sequence is given 
in Section 3. Using this algorithm we computed 10 7 terms; this led us to conjectured asymptotic 
formulae given in Section 4. We give a heuristic argument why these formulae may be true, but 
it seems likely they will be very hard to prove. We are able to rigorously extablish linear upper 
and lower bounds on the EKG sequence, namely < a(n) < 14n (Sections 5 and 6); the proofs 
use sieving ideas. Section 7 discusses experimental results concerning the cycle structure of the 
associated permutation of N. The final section discusses generalizations to other sorts of integer 
permutations resulting from greedy constructions with restrictions on gcd's of consecutive terms. 

2. The sequence is a permutation 

We begin with some general remarks about the sequence. 

For n > 2, let g = gcd{a(n — 1), a(n)}. For some prime p dividing a(n — 1), a(n) is the smallest 
multiple of p not yet seen (otherwise the smaller multiple of p would be a better candidate for 
a(n)). We call such primes p the controlling primes for a(n). There may be more than one, and 
their product divides g. 

For any prime p and n > 2, let B p {n) be the smallest multiple of p that is not in 
{a(l), . . . , a(n — 1)}. For example, the sequence {B2(n) : n > 2} begins 2, 4, 6, 8, 8, 8, 8, 10, 
14, .... Clearly 

B p (n) < B p (n + 1) < pn (1) 
for all p, n > 2. Then we have a(l) = 1, 

a(n) = mm{Bp(n) : p divides a(n — 1)} , (2) 
for n > 2, which provides an alternative definition of the sequence. 

Lemma 1. Let p be a prime > 2 that divides some term of the sequence. If p first divides a(n) 
then a{n) = qp where q is the smallest prime dividing a{n — 1), q is less than p, a{n + 1) = p, and 
either a(n) or a(n + 2) is equal to 2p. The new primes that divide the terms of the sequence appear 
in increasing order. 

Proof. Let a(n) be the first term divisible by p. The numbers pq where q is a prime dividing 
a(n — 1) are all candidates for a(n), and so a{n) = pq where q is the smallest such prime. Also p 
must be the smallest prime that has not appeared as a divisor of {a(l), • • • , a(n — 1)} (for if p' were 
a smaller such prime then p'q would be a better candidate for a{n)). In particular the primes that 
divide the terms of the sequence must appear in increasing order, and q < p. Then p is a candidate 
for a(n + 1), and is less than B q (n + 1) > B q (n) = pq, so a(n + 1) = p. Finally, either a(n) = 2p or 
else 2p is the winning candidate for a(n + 2). ■ 

Lemma 2. The primes that appear in the sequence occur in increasing order. 
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Proof. This follows from Lemma [T], since the first time p divides a term of the sequence the next 
term is p itself. ■ 

Lemma 3. If infinitely many multiples of a prime p appear in the sequence then all multiples of p 
appear. 

Proof. We argue by contradiction, and let kp be the first multiple of p that is missed. Choose 
no so that a(n) > kp for all n > hq. Since infinitely many multiples of p occur, there exists n> n$ 
with a(n) = Ip for some I. But now we must have a(n + 1) = kp, because gcd{a(n), kp} > p is 
allowed, and all smaller possible values which are ever going to appear in the sequence have already 
appeared. This is a contradiction. ■ 

Lemma 4. If all multiples of a prime p appear in the sequence then all positive integers appear. 

Proof. Again we argue by contradiction and let k > 2 be the first integer that is missed. Since 
infinitely many multiples of k occur among all the multiples of p, we get a contradiction just as in 
Lemma |3| Namely, there exists for n > hq a value a(n) = kip for some I, and gcd{a(n), k} > k is 
allowed, and all smaller possible values have already been used. Thus a(n + 1) = k, a contradic- 
tion. ■ 

Theorem 1. {a(n) : n > 0} is a permutation of the natural numbers. 

Proof. No number can appear twice, by construction, so it suffices to show that every number 
appears. Suppose only finitely many different primes divide the terms of the sequence. Then one 
of them would appear infinitely many times, and Lemmas |3| and [| would imply that all integers 
occur, which is a contradiction. 

Therefore infinitely many different primes p divide the terms of the sequence. Then by Lemma 
H infinitely many even numbers 2p occur, by Lemma ||| all even numbers occur and by Lemma [| all 
positive integers occur. ■ 

Remark. As will be discussed in Section 8, the principle of this proof generalizes to a wide variety 
of other integer sequences defined by restrictions on the gcd's of consecutive terms. 

3. Numerical investigations 

To compute the EKG sequence it is better not to use the original definition but to use dD and to 
store the current values of B p (n) for primes p. An efficient way to arrange the computation is to 
maintain four tables: 

hit(m) = if m has not yet appeared, otherwise 1; 

gap(m) = current value of B m (n) if m is a prime, otherwise m; 

small(m) = smallest prime factor of m; 

quot(m) = largest factor of m not divisible by small(m). 

Combining these tables in a C "struct" minimizes memory access. 

Suppose we wish to compute the sequence until a(n) reaches or exceeds N. The first step is 
to precompute small(m) and quot(m) for m < N. Since it is only necessary to consider primes 
< \/N, this takes about S p< ^N/p = 0(N log log N) steps. 
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In the main loop, let a(n) be the current value. Set k = a(n), B = N and repeat until k reaches 

1: 

p = small(A;) , 

B = m.in{B, gap(p)} , 

k = quot(fc) . 

Then we set a(n + 1) = B, hit(-B) = 1, and update gap(g) for primes q dividing B. 

We have not analyzed the complexity of the main loop in detail, but it also appears to take 
roughly 0(N log log N) steps, comparable to and not much greater than the number of steps needed 
for the precomputation part of the calculation. The program computed 10 7 terms of the sequence 
in less than a minute. For example a(10954982) = 11184814. 

4. A conjectured asymptotic formula 

The results from the experimental data suggest that whenever a prime p occurs in the sequence 
it is preceded by 2p and (consequently) followed by 3p. Although it is theoretically possible that 
some other multiple of p occurs before p, for example, we might have seen . . ., 3p, p, 2p, . . ., this 
does not happen in the first 10 terms. 

Conjecture 1. Whenever a prime p occurs in the sequence it is immediately preceded by 2p (and 
hence followed by 3p) . 

The numerical results also strongly suggest that the terms of the sequence fall close to three 
lines (see Figs. ||, ||). 

• if a(n) = m and m is neither a prime nor three times a prime, then a(n) ~ n; 

• if a(n) = p, p prime, then a(n) ~ n/2; 

• if a(n) = 3p, p prime, then a[n) ~ 3n/2. 

This was also observed by Ayres [Ayres 2001]. In fact, if we smooth the sequence by replacing 
every term a(n) = p or 3p, p prime > 2, by a(n) = 2p, the terms of the sequence lie close to a 
single line (see Fig. |5|). 
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Figure 5: The sequence smoothed by replacing a(n) = p or 3p, p prime > 2, by a(n) = 2p. 

A plausible but nonrigorous argument suggests a more precise conjecture. 

Conjecture 2. Let f(n) ~ g(n) mean that the ratio of the two sides approaches 1 as n — > oo. 
(1) If a(n) = m, m ^ p or 3p forp prime, then 

a(n) ~ n ( 1 + ) ; (3) 

\ 3 log n J 



(2) If a(n) = p, p prime, then 



a{n) ~ \n ( 1 + — — ) ; 
2 \ 31ogn/ 



(3) If a(n) = 3p, p prime, then 

a(n) ~ -n ( 1 + — — I • (5) 
W 2 V 31ogny v; 

To see why this conjecture might be true, consider a term a(n) =mof the smoothed sequence. 
Examination of Fig. ^ suggests that the smoothed sequence has hit all the numbers from 1 to 
m at least once (the numbers occur a little out of order, but never by much). However, we have 
smoothed away the numbers p and 3p that are < m, while picking up the primes p and 3p for 
2p < m. Therefore 

n^m-7T(m)-TT^—j+27T^—J, (6) 
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Figure 6: Values of c in (Q) for n near 10 7 . 



where tt(x) = number of primes < x. Then (||) follows at once from the asymptotic formula 
7r(x) ~ x/logx. Equations (f|) and ^) are based on (^) and the observations made at the beginning 
of this section. 

Although we are unable to prove this conjecture, it is an excellent fit to the data. 
If we try to write 

a(n) « n ( 1 + — L- + jr^) (?) (7) 
\ 31ogn (lognyj 

then the values of c do not appear to converge to a single value (see Fig. ||) , although c is very often 
close to 0.11. It seems conceivable that c might converge in distribution to a limiting distribution. 
Using two terms of the asymptotic expansion of ir(x) would give 



1 d 
a(n) ~ n I 1 + — h 



3 log n (log n) 2 



(?) 



(8) 



where d = 4/9 + (log3)/3 - log 2 = 0.1175. . . . 

Conjectures [l] and ^ predict that the fe-th prime p^ will occur in the pattern a(n) = 2pk, 
a{n + 1) = p^, a(n + 2) = 3pk, where 



2 Pk 



1 + 



31og(2p fe ) 



These conjectures may be hard to settle, because the permutation a(n) encodes an intricate in- 
teraction between additive and multiplicative properties of integers, which by the "greedy" property 
of the definition depends on all the earlier terms of the sequence. 
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In the next two sections we establish linear upper and lower bounds for the sequence, namely 

n < a(n) < 14n. 

260 - W - 

The numerical evidence supports the following conjectural bounds. 

Conjecture 3. The sequence a{n) satisfies a(n) > ^|n, with equality if and only if n = 28, and 
a{n) < Y~n, with equality if and only if n = 7. 

For large n the asymptotic lower and upper bounds would be n/2 and 3n/2, by Conjecture |2[ 

5. A linear upper bound 

In this section we show: 

Theorem 2. a(n) < 14n, for n > 1. 

Remark. The proof is by contradiction, and the basic idea of the proof uses the fact that the 
function values {a(k) : 1 < k < n} have a sandpile-like structure, in which each element can be 
built only at the top of a ladder of all smaller multiples of a controlling prime p dividing it. To 
reach a number at height exceeding 14n via a ladder of multiples of a controlling prime p, one 
repeatedly falls off this ladder while building it, at various smaller multiples kp of p of size between 
2n and 14n (see property (P4) below). To get back on this ladder one must use ladders of other 
controlling primes q which reach to some currently omitted multiple of p at the top of their ladder. 
The total number of elements in such ladders is shown to be large using a combinatorial sieve 
argument (compare [Hooley 1976], pp. 4-5); a contradiction results by showing that the sandpile 
contains more than n elements. 

We begin with a preliminary lemma. 
Lemma 5. If a{n) is divisible by a prime p, then p < n. If p ^ 2, p < n. 

Proof. The result is true if p = 2 or n < 3, so we may assume p > 3 and n > 4. We may also 
assume a(n) is the first term divisible by p. By Lemma |l[ a(n — 1) = pq where q is the controlling 
prime for a(n). Then pq = a(n) = B q {n) < q(n — 1) by ([!]). ■ 

Proof of Theorem ^|. We argue by contradiction, and let n be the smallest number such that 
a(n) > 14n. By direct verification we know n is large (in fact n > 10 7 , although we will not use 
that in the proof). Let p be the smallest controlling prime for a(n), say a(n) = Ip > 14n, so / > 15. 
Also a(n) = B p (n), so from (|l]) and Lemma [|, 17 < p < n. 
Let = \l/2~\ > 8, and consider the "window" 

w = {i'p, (i' + i) P , ...,(/- iM . 

Since p is a controlling prime for a(n), every element of W has already appeared in the sequence: 
a{t) G W implies t < n — 1. Call a(t) G W an entry point if a{t — 1) W, an exit point if 
a{t + 1) G" W. The following properties hold. 

(PI) The number of entry points is equal to the number of exit points. 
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(P2) At most one entry point a(t) has p as a controlling prime. This can only happen if a(t — 1) = 
(I' — l)p, a(t) = I'p, for some t < n — 1. 

(P3) At most one exit point a(t) has p as a controlling prime for a(t + 1). This happens just if 
a{n — 1) G W. Furthermore, if this is not the case then a(n — 1) = ip with i < and hence there 
is no entry point with p as a controlling prime. 

(P4) Let a(t) = ap G W, a £ [I', / — 1]. If a is a multiple of a prime q < 7, then a{t) is an exit 
point. (For a(t + 1) < 5,(t + 1) < 7t < 7(n - 1) < Z'p.) 

For a set of primes V, let D-p(a,b) denote the number of integers a < i < b such that i is a 
multiple of some element of V . Let 

e = D {2A5jT} (i',i-i). 

By (P4), the number of exit points is at least 9, while if there exists an entry point controlled by 
p then there are in fact at least 9 + 1 exit points. By (PI) we conclude that the number of entry 
points not controlled by p is at least 9. 

Let S = {qi, (/2, • • •} (with q\ < q 2 < ■ ■ •) denote the set of controlling primes for entry points 
not controlled by p. Note that 2,3, 5,7, p are not in S (by the same argument as in (P4)), and 
\S\ < ir(l - 1) - 4. Then 







< 
< 

< 
< 



£> {2 ,3,5,7}(^-l) 

number of entry points not controlled by p 
D s (l',l-1) 

E r«-r- 

qes 



(i-i'-i)yl + \s\. 



Setting § = ^ qeS h we have 



> 



tt(Z - 1) +4 



i 



(9) 



(10) 



The right-hand side of (|T^) is a function of the single variable /, and is > 2/9 for all / > 15, with 
equality if and only if / is 20 or 21. (This is easily verified by computer for small I, say / < 1000, 
and analytically for larger /.) In other words, 



x ^ 1 2 

4^% " 9 



i>l 



fc-1 



Define k by 

i=l % y i=l q l 

and let S' = {qi, . . . , qk} Q S. Note that q^ > 17, since 1/11 + 1/13 + 1/17 > 2/9, but no proper 
subset of {1/11, 1/13, 1/17} has this property. 

Every multiple of any element q G S that is < I'p must have already occurred in the sequence, 
by definition. We obtain a contradiction by showing that there are more than n different multiples 
of elements of S that are < I'p. By inclusion-exclusion we have 

k 



Ds(l,l , p)>D S '(l,l , p)>J2 



Vp 



E 

l<i<j'<fc 



I'p 



(12) 
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To bound the first term in (12), observe that for qi £ S there is a multiple of qip in W, so 



14n Ip 



E 

i=i 



7p 



> 



> 



8Z'p 
9^7 



> ° V — 



To bound the second term in (12) we use 



E i 



i<i<j<k 



^fe-i , \ 2 , fe-i 



< 



2 \ " q 

\i=i H 



2Qk f-f ft 



Then, since > 17, we have 



fc-l 



2 1 /2 



9 9 2 V 9 



1 2 
34 ' 9 



229 1603 

> In = n > n , 

1377 1377 



which is the desired contradiction. 



6. A linear lower bound 

In this section we show: 

Theorem 3. a(n) > [^Bo^l; for n > 1. 

Remark. The proof is a modification of that of Theorem |l|. It aims to show that if some number 
less than ra/260 is missed in {a(k) : 1 < k < n} then there are at least n/65 numbers in this 
set that are even numbers, and Lemma [?] below provides the mechanism to get a contradiction. 
The method of Theorem |l| seems inherently weaker when used for a lower bound, so we have not 
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attempted to streamline the proof. It would certainly be possible to reduce the constant 260, but 
not to anything close to 14. 

We begin with three lemmas. 
Lemma 6. For a prime p, if a[n) = kp for some k then a(j) = k for some j < n + 1. 

Proof. We argue by induction on k. The result is true for k = 1 since a(l) = 1. 

Let q be a controlling prime for a(n), so q\a{n — 1) and q\a{n) = kp. Case (i): q ^ p. Then 
q\k, say k = mq, a(n) = mqp. Hence all multiples iq with i < mp have already appeared, and in 
particular a(j) = mq = k for some j < n. 

Case (ii): q = p. All multiples ip with i < k have already occurred. By the induction hypothesis, 
i has occurred for all i < k. If k has occurred then a(j) = k with j < n, and otherwise a(n+l) = k. m 

Lemma 7. // at least 4k even numbers occur in {a(l), . . . , a(n)} then all numbers {1, . . . , k} occur 
in {a(l), . . . , a(n + 1)}. 

Proof. In view of Lemma || (taking p = 2) it is enough to show that {2,4, ... ,2k} are in 
(a(l), . . . ,a(n)}. 

Suppose not, and let 2m be the largest even number < 2k not in {a(l), . . . ,a(n)}. Every even 
number a{i) > 2m with i < n will be followed by a{i + l) < 2k (since 2m is always available). But in 
{a(l), . . . , a(n)} we have at least 4k even numbers, and so at least 2k even numbers > 2k. Therefore 
in {a(l), . . . , a(n + 1)} we see all the numbers from 1 to 2k, including 2m, a contradiction. ■ 

Lemma 8. If a(n) = p, a prime, then all numbers {1, ... ,p — 1} occur in {a(l), a{n — 1)}. 

Proof. By Lemma |l|, a{n — 1) = qp, q prime, q < p; so q, 2q, . . . , (p — l)q have already appeared 
in {a(l), . . . , a(n — 2)}. The result now follows by Lemma ||. ■ 

Proof of Theorem By direct verification we may assume n > 260 2 . Let m = [n/260], and 
suppose, seeking a contradiction, that a(n) < m. Note that the lower bound on n implies that 
n/m > 259. 

We will show that at least 4m even numbers have occurred in {a(l), . . . , a(n — 2)}, which gives 
a contradiction by Lemma |?|. No primes greater than m can occur in this interval, or we get a 
contradiction by Lemma [8|. 

Some number > n must occur among {a(l), . . . ,a(n — 1)}, since there are n — 1 numbers and 
a(n) < m is missing. Suppose a(j) > n, with controlling prime p < m, say a(j) = Ip, and j < n—1. 
Since Ip > n, I > 260. Let R < m be the smallest missing number among {a(l), . . . ,a(n — 1)}. 
Then / < R, for if / > R then we have seen Rp at time j — 1, and by Lemma |6| we have seen R 
by time j, a contradiction. Therefore I < m and so p > n/7 > 259. Both / and p are in the range 
[260, m). 

We consider the "window" 

+ p,...,lp\ , 

and define "entry" and "exit" point as in the proof of Theorem |2[ 

There are at least [3Z / 4J multiples of p in W, and at least [3Z / 8j — 1 of them are even. Any such 
even multiple of p is an exit point (from Lemma |?]) . There must therefore be at least [31 /8j — 2 
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entry points which are not controlled by p. Let S = {q\, q2, ■ ■ •} (with q\ < q<i < ■ ■ ■ < I) be the set 
of controlling primes for these entry points. Then 



As in (|), (0), we get 



E^ 

q% 



lj > [3l/8\ - 2. 
LfJ -2-vr(Z) 



1 LiJ 



(13) 



which is > 43/214 for I > 260. We define k by 



k-l 

E7 



i=l 



43 

< < 

qi 214 " 



k 

E^- 



(14) 



i=l 



Every multiple of any element q £ S that is < [Zp/4J must have already occurred in {a(l), . . . , a(n- 
2)}. Of these, at least [lp/8q\ are even. Therefore the number of distinct even numbers in the 
range \}p/&\ that have occurred is at least 



E 

i>l 



X 
2q~ 



E 

l<i<j<fc 



X 



where X = Ip/A. Proceeding as in the proof of Theorem |2|, and using q\ > 2, we find that this is 
at least 4.024m. This exceeds Am and provides the desired contradiction. ■ 



7. Cycle structure 



Since the sequence is a permutation of N, it is also natural to investigate the cycle structure. 
The experimental evidence suggests that there are infinitely many finite cycles and infinitely many 
infinite cycles. It seems very likely to be hard to prove either of these observations, or even to prove 
that there is at least one infinite cycle. 

The first few finite cycles start at the points 

1, 2, 3, 8, 40, 64, 121, 149, 359, 2879, 5563, 28571, 251677, 

and have lengths 

1,1,6,1,1,1,2,12,11,25,8,22,11 
respectively. There are also a large number of apparently infinite cycles, of which the first two are 

• • • 229310, 117833, . . . , 22, 27, 26, 28, 13, 14, 7, 12, 18, 20, 11, 15, 21, ... , 636551, 652766, . . . 

and 

. . . , 502008, 257519, . . . , 248, 253, 131, 138, 73, 82, 129, 201, 212, . . . , 645906, 662330, . . . 

with minimal representatives 7 and 73 respectively. The first fifteen of these apparently distinct 
cycles have not coalesced in the first 700000 terms of the sequence. However, although it seems 
unlikely, it is theoretically possible that they could coalesce at some later point. It would be nice 
to know more! 
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8. Generalizations 



The EKG sequence can be generalized in various ways, while retaining the basic construction of a 
greedy sequence with a condition on gcd's of consecutive terms. For fixed M > 2, let b(n) = n for 
1 < n < M, and for n > M + 1 let bin) be the smallest natural number not already in the sequence 
with the property that gcd{6(n — 1), 6(n)} > M. The proof of Theorem || easily extends to show 
that (b(n) : n > 1) is also a permutation of N. For the cases M = 3,4,5 see sequences A064417 , 
A064418j |A064419| in [Sloane 2001]. 
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Figure 7: Plot of 6(1) to 6(1000) for the case M = 3. 

The first 1000 terms for the case M = 3 are shown in Figure ^. This sequence appears to 
behave in a similar way to the EKG sequence: the spikes in the sequence are associated with the 
primes, occurring in the order 3p,p, 2p,4p. All points of the sequence seem to lie near lines of slope 
1/3,2/3,1 and 4/3. 

More generally, we can allow an arbitrary finite prefix at the beginning of such sequences. 
Consider a sequence c(n) with the generating rule that gcd{c(n — l),c(n)} > M, for n > N, but 
with a finite prefix c(n) = a n : 1 < n < N, such that all terms a n are distinct natural numbers, 
and that the prefix includes values a(n) = k for of 1 < k < M. With a little work the proof of 
Theorem || can be modified to apply to these sequences as well to show they are permutations. 

If we use the greedy gcd rule with M = 2, starting with an arbitrary finite prefix, one obtains 
an infinite number of different EKG-like permutations. Many of these sequences, perhaps all, 
eventually coalesce with the original EKG sequence. In any case, the qualitative properties of 
the resulting permutations appear similar to the EKG sequence: the primes provably appear in 
consecutive order (excluding those in the prefix terms); the general plot of the permutations remains 
similar to Figure ||. 

It appears that linear upper and lower bounds hold for all the permutations c(n) of these types. 
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Very likely such bounds can be established rigorously in any particular case using methods similar 
to those used in proving Theorems [2] and ||. 
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